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1.  INTRODUCTION 


The  Radiation  Transport  Committee  was  asked  to  consider  th r 
following: 

1.  I'  sonant  scattering 

2.  Excitation  transfer 

3.  Self  -reversal 

4.  Atmospheric  transmission 

5.  Cloud  effects 

6.  Albedo 


7.  Factors  that  determine  time  dependence  of  signal 

8.  Solar  scintillation 


The  study  of  resonant  scattering  and  excitation  transfer  for 
the  Lyman- Birge-Hopfield  series  of  nitrogen  and  for  the  first  posi 
tive  system  of  nitrogen  has  been  treated  by  Ed  Zipf,  whose  report 
is  included  in  Appendix  A.  Self-reversal,  which  occurs  when  an 


optically  thick  medium  is  under  consideration,  is  not  important  to 
this  discussion  because  th* limiting  range  for  detection  of  an  existing 
system  is  one  in  which  the  X-ray  excitation  is  so  mild  that  the  fluor¬ 
escent  layer  behaves  an  an  optically  thin  medium. 

A  general  discussion  of  atmospheric  effects  pertinent  to  the 
VELA  program  is  presented  below.  The  material  in  Section  2  on  atmos 
pheric  refraction,  absorption  and  scattering  in  the  atmosphere,  sky 
backgrounds,  noise  in  the  sfty,  and  reflectivity  of  the  earth's  surface 
is  taken  from  Ref.  1.  The  material  on  typical  probabilities  of  signal 
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modifications  by  clouds  and  overcasts  is  taken  from  Ref.  2.  Work  on 

th*  Carth'S  albed°'  a£  satellites,  and  the  effects  of 

this  on  background  and  time  delay  of  signals  was  summarised  by 

L-  F.  Drummeter  and  incorporated  in  the  Satellite  Committee  Report. 
The  material  presented  in  Section  2  on  the  delay  of  light  pulses  in 
Overcasts  and  clouds  has  been  abstracted  from  a  series  of  reports 

prepared  at  E.  H.  Plesset  Associates  under  AFTAC  Contract  No. 

AF  ?3(657).  13599.  The  application  of  these  delaying  effects  to  the 
specific  VELA  problem  is  illustrated  in  Appendix  B. 

The  work  on  scintillation  by  G.  Milne  is  presented  in 
Appendix  C.  The  scintillation  of  an  extended  source,  calculated  at 

E.  H.  Plesset  Associates  under  AFTAC  Contract  No.  AF  33(657)- 1 3595, 

U  PreS°'*ted  APP“**  D  and  summarized  in  Section  2  below.  Exper-  ’ 

imental  work  on  the  scintillation  of  the  sun,  done  under  the  same 

AFTAC  contract  at  E.  H.  Plesset  ARc^;,f 

Plesset  Associates,  is  reported  in  Appendix  E. 


0 


0 


2.  ATMOSPHERIC  EFFECTS 

Much  of  the  material  on  refraction  and  absorption  is  taken  from 

3 

Allen.  An  excellent  review  of  atmospheric  refraction  is  given  by 
Mahan.  ^ 

2.  1  Atmospheric  Refraction 

The  earth  is  a  sphere  of  radius,  R  =  6.  37  x  103  km.  The 
distance  D  t£>  the  geometrical  horizon  for  an  observer  at  height  h  («R) 
l8D  ^  2Rh  °r  D<km)~3.8  /h(meters)  or  D(nautical  miles)~l.i3 
\/h(ft).  A  target  at  h'  can  just  be  seen  over  the  horizon  by  an  observer  at 
h  and  at  a  distance  D*  =  >/2Rh  +  V2Rh\  Refractive  effects  can  modify 
these  results  by  approximately  +  20%  in  D. 

Values  of  (n  -  1)  x  10^  for  various  A  and  for  p  =  760  mm 
t  =  0°C  and  f( water  vapor  pressure)  =  4  mm  Hg  are  given  by  Table  I. 


Table  1  --  Index  of  Air  as  a  Function  of  Wavelet 


th 


- g - P-20  0.  24  0.  28  0.  32  0.  36  0.  aq  q.  5q  qj6o  0  10  QeQ 


(n-l)xlO  340.0  321.2  310.9  304.9  300.9  298.2  294.1  292.0  290.7 


290.0 


0.  90  1.  0  1.  2  1.  4  1. 6  1.  8  2.  0  3.  0  4.  0 


oo  optical 


2M 

(n-l)xlO  289.4  289.0  2886  2883  2881  2880  237.9  287.7  2S7.6  287.T"~ 

I 

The  combination  of  the  altitude  gradient  in  refractive  index  and  the 
earth's  curvature  result  in  a  deviation  of  the  light  from  stars  and  a 

limitation  on  the  air  mass  ’  e tween  star  and  ground  observer  as  shown 

» 


Trap  r.euith 
diauaiue 


Apparent 
altitude 
h  rn  00  - ; 


Table  II  —  Refraction  and  Air  Mass 

Refraction  ami  air  maw.  The  table  -l-nv,  tho  relation  between  refraction  ar  ,i  A.r  .. 
?.cn»tb  distance  s;  it  also  gives  tlic  air  mass  which  fo-  smnti  •  , 

Tho  values  arc  tabulated  for  p  ~  TOO  mm-H-  and  t ..  ‘in  *c  ‘‘r  ^  *°  6'C; 

i£ZfrUh,ly  both  tho  L  ~  c-X°;xuoc 


Crude  calculations  of  altitude  dependence  of  pressure  to  about  120  k.- 
ntay  be  made  using  the  height  of  the  homogeneous  atmosphere  H  as  a  value  fo 
scale  height  ?h  ■=  ?c‘h/H  where  H  is  related  to  ground  temperature  as  shown 


in  Table  in. 


Table  HI  — 


Height  of  the  Homogenoous  Atmosphere 


Ground  Temp.  °C  -30  .15  p  tl5  t 

H  in  km  j  7. 113\  7. 552  7.991  8.  43  8.  So9 


Reasonable  calculations  ofjaltitude  dependence  of  temperature,  a  uo 
to  the  tropopause  may  be  Aade  as  follows:  \  =  t  .  (t  -  .  j  y_/\  ' 

The  variation  of  tropopausi  height,  h,  and  •  «emp0eratu°ro  t  with'’ 
changing  latitude  is  shown  la  Table  IV".  C 
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A  '  \ 


•.bt’  .n.v  21* 

S!>-  24'  40* 
SS*  IS'  24' 
S7  14'  24* 
S«3J  11' 43* 
SO*  0'51* 


S3*  "  23 
S2*  6' 33* 
SI*  5' 52* 
SO*  5'  IS* 


Apparent 

altitude 


Truo  zenith 
distance 

c 


Air 
m.m 
[3.  01 


Beer 


so  y  is* 

75*  3'  34' 
70*  S’  3S* 
CO’  2'  4* 
CO*  1'  41* 

50*  r  io* 

40*  0'  40* 
30*  0'  34' 
20*  0'  21' 
10*  0'  10* 
0*0'  0* 


5-700 

3-SC4 

2-024 

2-30C 

2-000 

l-e5C 

1-305 


O.f 


•  'Od  I 


1-000 

1-553 

1-304 

1-154 

1-00-i 

1-015 

1-000 


1-0C4 

1-015 

1-000 


Table  IV  --  Height  of  the  Tropopausc 


Variation  of  Meteorological  Quantities  with  Latitude 

The  table  averages  (ho  Northern  arnl  Southern  hemispheres  which  differ  in  dev 
as  a  result  of  the  different  land  distributions 


Loti* 

tudo 

Mean  air 
temperature 
(.ca-lcvcl) 

Seasonal 
tempera  turo 
rango 
(on  land) 

•c 

*C 

0* 

27 

1 

10* 

20 

3 

20* 

24 

C 

30* 

20 

'  0 

40* 

13 

13 

SO* 

•  6 

17 

CO* 

—  *> 

21 

70* 

-10 

2C 

SO* 

-IS 

20 

00* 

-25 

Total  Water..  •  Tro?c, 

pressure  V£“I>our  - — — 

(sea-level)  temperature  )  hc,XM 


mm -Kg 

75S 

750 

7CI 

7G3 

7C1 


mm -Kg 


2.  2  Absorption  in  the  Atmosphere 

The  specular  transmission  through  a  homogeneous  atmosphere  Ti  of  mono- 

*  t  *  ' 

chromatic  light  of  wavelength  A  over  equivalent  path  length  L  in  cm  at  ST? 
is  given  by 

t\  *«p-(K)k+rx)i.  '  U) 

w.nere  is  the  absorption  coefficient  and  (T^  is  the  scattering  coefficient. 
Sma.il  wavelength  intervals  in  the  ultraviolet  and  visible  may  be  treated  by 
r*q.  (i)  and  values  of  are  given  in  Figure  1.  Appropriate  values  of  L  may 
be  derived  from  Figure  2.5  In  cases  of  N2,  02,  03  and  H20  the  composition 
of  unit  atmosphere  in  cm  thicknesses  at  STP  are  624,  600,  167,  600,  0.  3  and 
from  300  to  5000  respectively. 

In  the  infrared  the  rotational  structure  of  the  absorption  spectra  c 
atmospheric  molecules  is,  in  general,  not  resolved  and  the  assumotior.  tha 

■  i 

is  nearly  constant  over  a  stnall  wavelength  interval  is  not  valid.*  Ur.de 
these  conditions  the  transmission  is  described  by  Eq.  (2). 

Ta  =  f(bA<0  exp(-<r.  L)  j 
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Figure  1  (from  ref.  3) 
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s*c  i(b^ ^  )  is  given  m  1  able  V  and  values  for  log  b^  arc  shown 
s  3  and  4.  Here  J!  corresponds  to  the  values  of  b\  for  unit  air 


Table  V 


Band  Absorption  of  Atmospheric  Cases 


llrtmU  matte  vip  of  iliac  reto  lino*  will  not  obey  tlio  Lacnkeri  low  of  exponent  L 
tion  and  the  absorption  coefficient  must  be  rcpi«C'*d  by  it*  analogues  i»A  sue!; 
transmission  is /(/»*/)  whore  l  is  the  layer  thickness  in  soino  appropriate  tint 
function  /we  adopt  tho  relation  given  below: 


Except  for  the  Chappuis,  the  Huggins,  and  the  Hartley  bands  of  CL  zl 
little  signiiicant  absorption  in  the  atmosphere  between  0.  3  ji  and  0.  7 
Table  VII  ). 


log  (V) 

/(Ml 

0d>" 
+  0-2 

T  0’** 

+  0-C 
+  03 

+  10 

' 

0-000 

0-412 

0-332 

0-203 

0-137 

0-124 

,s  ftnioimc 


j  i  is 


IA  (ift  cm-*) 

*  v 

Injru-rcd  band  absorption  oj  atmospheric  gases 


Figure  4  (from  ref.  3) 


—The  efficiency  factor,  £t  lor  extinction  by  email  iphcrea  ae  a  function  of  the 
parameter,  ay  and  the  refractive  Index,  m. 


Figure  5  (from  ref.  7 ) 
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Figure  6  (from  re 
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Figure  3 
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Tabic  VII  --  Continual  Atmospheric  Absorption 


Molecular 

•cattcrinp; 


WuUtr- 
Vf pour 


Oiono 


fairly  clear 
comlitir.r.u 


per  cm 


--i*™  nr;b“ 


per  3  tr.m 
ftt  JJ.T.1*. 


7-20 
4-CS 
3- 1(J 

j.jj 
1-006 
H04 
0-00$ 
0-702 
0-652 
...  0-441 


per 

atmosphere 

0-24 

•21 

•19 

■17 

•157" 


0-356 
0-219 
0-142 
0-005 
0-0675 
0*04  SS 
0-0301 
0-0211 
0-0131 


0-00S6 
0-0041 
0-0022 
0-C013 
0-01 80 
0-0  53 
0-0*  10 
0-0‘33 


•0010 

•0010 

•0010 

•0010 

0-0010 


The  atmosphere  scatters  sunlight,  moonlight  and  starlight 
for  wavelengths  in  excess  of  0.  3 ^microns  (the  ozone  cutoff)  there  is 
varying  brightness  over  the  whole  sky.  In  addition  to  this  scattered 
is  a  background  Observable  at  night  which  is  due  to  other  sources  as 


Zodiacal 


Galactic 


Luminescence  of  night  sky 
Scattering  from  the  above 


Plus  starlight  direct  and  scattered 


•  * 


M. '  ~ 

f  x  I  |  j  <j 

.  1 


‘  P 


.  .  • 

■  ■> 


Rayleigh  scatto: 
ons  :'or  Rayleigh  scarce 


i*mg  causes  polarization  of 
!?iag  of  unpolarizcd  light  r: 


"  "  'Tott  cos  y  'vhcr«  J-  and  11  describe  the  st; 
no  of  scattering.  Hence  the  light  scattered  into  u: 
om  unit  volume  containing  m  scattering  centers  j 


or  the  molecular  scattering  in  the  atmosphere 
nd,  with  the  reservations  pointed  out  by  van  de 

X0  "isl - U'4  cos2  ®)  (1.  07  x  10~8/\  "4*  05)  /ne 


station  is  illustrated  in  Figure  9 


SUyWjh 

iaeidoat  rtduUon  i 
X  ?l*ao  of  dritrla; 


~  — v.o  wee  atmosphere  in  general  produce 
•-g  *-r.d  some  polarisation]  The  phase  functions  for  e, 
e  base  and  for  clouds  arejshown  in  Figure  10  a,  b.  a: 
oral: ration,  the  sky  brightness  above  the  tropopause 
■SiS  o!  “V-6  ^yleiglt'  scattering  phase  function  inereas 
:s  sometimes  holds  above'an  inversion,  and  infrce.ue; 
dr  the  conditions  in  which  the  single  seattert*.  _ _ 


Normalized  scattering  phase  functions 
atmospheric  conditions.  (Note:  to  find 
independent  values,  take  average  of  tw 
rather  than  sum). 


HAZE  MODEL  C 
X  0.45M 


-Vi.,'-'’ 


MODEL  V. 


+  0.4:8 


CLOUD  MODEL 
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acceptable  (total  path  length  in  atmosphere  less  than  cr  "  *)  spectral 

sky  brightness  values  are  given  by  I(  0,  z,A)  =  1^^  H(  see  z)  f(g) 

watt  cm  sr  ji  ,  where  Iq^j  is  incident  sunlight  in  watts  cm'2, 

2  iS  anSl°  £r°m  Zenith  to  Um  01  0  is  angle  through  which  sunlight 

is  scattered  toward  detector. 

The  transmission  or  reflection  coefficients  of  uniform  overcasts 
and  the  reflection  characteristics  of  large  plane  surfaces  of  clouds  can 
be  treated  with  considerable  accuracy  through  the  use  of  the  following 

equations.  The  transmission,  T,  of  the  overcast  for  light  incident  at 
an  angle  0  to  the  overcast  normal  is 

T  =  1  +  <cos  0  ~  »5)(1  -  e~P/A  cos  °) 


where 

D  is  the  thickness  of  the  overcast 
A  is  the  $ransp’ort  mean  free  path, 

A  =  V(17V* 

Ag  is  the  scattering  mean  free  path  in  the  overcast, 

Fo  is  the  average  cosine  of  the  scattering  angle. 

Hence,  for  light  incid  rt  at  an  angle  0  =  60°  the  transmission  is 

t6o°  =  i/i+  (pM )  (3) 

4L  T>  •  S: 

\ 

The  normalised  distribution  in  angle  of  the  transmitted  light  is  made 
up  of  two  parts.  One  part  je  from  the  light  which  was  scattered  within 
the  overcast.  The  other  isjfrom  the  non-scattered  light,  which  is 

c.ular  component  is  proportional  to  e“^/^ scos  ® 


E,  H.  PLESSET  ASSOCIATES,  INC. 


The  scattering  in  an  overcast  is  very  forward,  i.  e.  , 


A  »  A_  . 


Therefore,  even  for  TAfl0  =  .  5  (a  rather  thin  overcast),  the  factor 


D/A  is  quite  large;  hence, the  unscattered  component  may  be  ignored, 
s 


In  this  case  the  distribution  in  angle  of  the  transmitted  light  is 


0)  =  { 3/7 rr) ( 1+  2  cos  f) 


where 


0  is  the  angle  to  the  normal  ol  the  overcast 

f(^)  is  the  normalized  distribution,  l/steradian,  i.  e.  , 


fXj  2 

2  iri  0).  cos  sin  dj&  =  1 


Also,  the  amount  of  light  reflected  is  given  by  R  =  1  -  T,  or 


(D/A)  -  (los  0  -  ,  5)(i  -  e~D/A  COS  °) 

=  - p-r+ron 


And  this  light  is  distributed  in  angle  according  to  0). 


The  transmission  of  light  as  a  function  of  time  was  calculated 
using  diffusion  theory. 


Shown  in  Figure  His  the  probability  of  transmission  of  the  scat¬ 
tered  light,  P,  by  an  overcast,  as  a  function  of  time  for  various  over¬ 
cast  thicknesses.  The  source  is  an  instantaneous  pulse  of  radiation 


incident  at  an  angle  of  60°  to  the  normal  at  t  =  0.  The  probability  of 


transmission  of  the  scattered  light  as  a  function  of  time  is  normalized 
bi  that  the  integral,  P,  over  all  time  is  as  follows: 


/*“  -D/A  cos  0 

-  J  P  dts+  e  ' 
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igure 


Probability  of  Transmission  of  an  Overcast  Vs.  Ti 
for  Various  Overcast  Thicknesses 


an  instantaneous  pulse  incident  at  60°  to  the  normal 


Maximum  Probability 
One -Ha  if  Maximum 
One -Tenth  Maximum 


D  =  Thickness  of  Overcast 
A  .  =  Trancpvit  Mean  Free 


Time 


t:k  r, 


As  can  be  seen  in  the  figure,  the  shape  of  the  intensity  vs.  time 
distribution  is  the  same  for  clouds  where  D  » . 


The  time  of  the  maximum,  t_n,  in  the  probability  of  transmission 
is  closely  given  by 


c  is  the  speed  of  light 

A  line  through  the  points  corresponding  to  the  probability  at  the  time 
of  the  maximum  is  shown  in  the  figure. 

An  equation" for  D/A  m»y  be  obtained  from  Equation  (3),  i.  e.  , 


And  the  time  of  the  maximum  is  as  follows 


Hence,  the  time  of  the  m<odmum  is  a  factor  greater  than  one  multiplied  by 
the  time  light  takes  to  perietrate  straight  through  the  overcast. 


The  times  at  which  the  probability  reaches  a  given  proportion 

-3/2 

of  the  maximum  value  are  found  to  vary  as  t  '  for  rather  thick 
overcasts,  i.  e.  ,  D  »A  .  Lines  corresponding  to  probabilities  of 
i/10,  1/2  and  equal  to  the  maximum  value  are  shown  in  the  figure  and 
have  slope  of  roughly  -3/2.  The  times  at  i/lO  maximum  value  may  be 
expressed  in  terms  of  '•he  times  at  the  maximum  as 


and  the  times  at  1/2  maximum  value  are 


mtrared,  a.  wavelengths  greater  than  about  3  r 

good  approximation,  ihiclt  clouds, radiate  as  blach  bodies.  Th- 

^  TaMe  "  **  «"  —mission  *  hazes  im?rc 

length  increases  and  thus  th*  ♦■Vea«  i 

a  thus  the  thermal  emissivity  will  decrease 

wavelength.  At  10  microns  mnct  v 

n«  most  hazes  are  essentially  transpa 

tms  wavelength,  clouds  are  black  In  -wm-*. 

oxacit.  In  addition  to  the  emission 

droplets  the  atmosphere  emits  in  those  wavelength  regions  whe 

“  &tm°Sphere  abSOrK  3  -  «  ‘*ow  that  such  emissi 

“  "  regl0aS  batWa°a  th°  atmospheric  windows.  », 

7  £ly  “d  te  **  “  ‘  *  wavelength  is  presented  in 

3  respectively.  .  • 


cation  of  Figures  12and  13 


-  Scattered  sunlight.  V  =  94  km,  0  =  10°, 

# 

-  Scattered  sunlight,  V  =  94  km,  0  =  40°, 
■  50%  transmitting  cloud  or  fog  0  a  30°,  s 

Mattered  sunlight,  V  -  ?4  km,  0  =  90°, 

Scattered  sunlight,  V  =  94  km,  0  =  120°, 

Blackbody  at  283°  k 
€ 

Emission  of  water  vapor  and  CO 

2t 

Bright  Aurora 

Ha?e  Radiance  plus  scatter  of  earth  flux 
V  *  50  km,  z  a  60° 

As  "I",  but  V.  R.  =  04  w 


A  -  Blackbody  at  283°  k.‘ 

B  -  Emission  of  water  vapor 
C  .«*  Bright  Aurora 
Dr  Air  glow 

E  -  Haze  radiance  plus  scatt< 
E  -  As  "E”,  but  V  a  94  km 

f- 

G  -  Full  moon,  V  a  94  km,  0 
H  -  Full  moon,  V  =  94  km,  0 


F ull  Moon,  V 


mum  scattered  city  lights 

moon,  V  a  04  km.  &  =  in9 


0.  is  the  angle  through  which  light 
the.  viewing  path 
a  is  the  "  *nith  angle  of  the 


.owing  path 


2.  5  The  Effects  of  Clouds  and  Overcast 

The  effects  of  clouds  end  overcast  -ire  twofold :  First,  t hoy 
may  reduce  the  magnitude  of  the  signal  by  scattering  light  from  the  overcast 
tack  into  space,  and,  second,  they  may  distort  the  signal  which  would  have 
been  received  in  the  absence  of  tpe  crouds. 

,* 

Data  taken  by  Eastman  Kodak  Company  from  1037  to  cate  indicate 
that  in  Rochester,  New  York  the  probability  that  illumination  at  the  surface 
of  the  earth  will  be  cut  by  clouds  to  5  10  of  clear  weather  illumination  is  loss 
than  2.5  x  10  and  that  it  will  bo  cut  to  10  &  is  about  10  , 

Probability  that  there  will  be  some  signal  reduction  and/or 

distortion  due  to  clouds  and  overcast  have  bean  calculates  from  United 

States  Department  of  Commerce,  Weather  Bureau  Paper  No.  12,  “Sunshine 

and  Cloudiness  at  Selected  Stations,  1951",  and  from  A?  Weather  Service, 

H 

Directorate  of  Climatology,  Data  Control  Division's  reports,  "Uniform 

i* 

Summary  of  Surface  Weather  Observations."  The  results  are  given  in  Table  VIII 


I  ~27' 

i 

I 


Total  probability  of  some  signal  nicdificatior 
or  clouds  (in  per  cent) 


p 


Z.  51  Cloud  Characteristics 


The  following  information  about  cloud  characteristics  was 


taken  from  the  ''Handbook  of  Meteorology  by  Berry,  Bollay 


Beers*  McGraw-Hill  (1945) 


Clouds  are  divided  into  four  families  with  reference  to  height 


above  the  ground  level  in  the  region, 


High  clouds  (mean  lower  level  6  km) 

Middle  clouds  {mean  upper  level  6  km,  mean  lower 
level  2  km) 

Low  clouds  (mean  upper  level  2  km,  mean  lower 
level  close  to  the  ground) 

Clouds  with  vertical  development  (mean  upper  level 

that  of  cir-rus,  mean  lower  level  . 


Family  B 


Family  C 


Family  D 


At  nearly  all  levels,  clouds  may  appear  under  the  following  forms: 

a.  Isolated  heap  clouds  with  vertical  development  during  thei 
formation,  and  a  spreading  out  when  they  are  dissolving. 

b.  Sheet  clouds,  which  are  divided  into  filaments,  scales,  or 


rounded  masses,  and  which  are  often  stable  or  in  process  of 


disink  gration 


c.  More  or  less  continuous  cloud  sheets,  often  in  process  of 


formation  or  growth 


The  following  table  indicates  the  cloud  genera,  with  the  observed  form 


■7* :  - 
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Family- 


Genus 


Cirrus 

Cirrocumulus 
Cirro stratus 


Altocumulus 

Altostratus 


Stratocumuius 

Stratus 

Njmbostratus 


Cumulus 

Cumulonimbus 


■ 

mm  .  ay 
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r*  WIT 
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■  Under  ideal  conditions  and  with  a  perfect  -adiation  detector  the 
noise  in  a  signal  or  the  noise  from  a  background  can  be  calculated  on  the 
.  basis  of  statistical  fluctuations  in  the  number  of  photons  reaching  the  detector. 
This  situation  almost  never  occurs  in  atmospheric  work.  Noise  in  me  signal 
arises  from  scintillation  and  seeing  problems,  and  noise  in  the  background 
lavel  .«  generated  by  space  and  time  variations  in  the  scattering  or  emitting 
media  responsible  for  background  signals.  Particular  interest  in  the  back- 
ground  noise  seen  by  systems  which  scan  the  sky  has  resulted  in  studies  of 
sky  noise  as  a  function  of  scanning  rate.  Resulting  data  are  usually  summer- 

_  1Md  ia  4  ta°wn  *»  *  Wiener  Spectrum1.2  Pour  typical  Wiener 

*  Spectra  for  a  circular  scan  are  shown  in  Figure  14. 

Turbulence  in  the  atmosphere,  particularly -in  the  tropopause  causes 
spatial  and  time  variation,  in  refractive  index.  One  may  thuds  of  this  effect 

"  pr0dUCi"2  refr4CtlVe  C0U*  tensions  of  the  class  of  three  to  eight 

inches.  Cells  deviate  the  Ught  passing  through  them  causing  variations  in 

the  position  of  the  stfciar  image  in  a  telescope  (seeing)  and  also  fluctuating 

shadows  on  the  aperture  of  the  telescope  (scintillation)..  The  lifetime  of  the 

ceils  *s  oi  the  order  of  0.  OS  seconds  and  they  may  move  with  the  winds  aloft 
(velocity  up  to  120  mph). 

For  observatories  on  the  earth's  surface  astrenemicai  resolution 
ts  limited  by  the  seeing  described  above  to  approximately  3  x  10'4  radians. 

&  *er-er41,  d4ytime  re£0luti0“  does  not  exceed  1  or  2  X  10'5  radians  for  a 
vertical  path.  Approximate  calculations  can  be  made  for  camera  systems 
looking  down  through., he  atmosphere  using  the  conation  A  x  .  cL .  .  h.  where 
A  x  =  minimum  resolvable  distance  on  ground,  measured  perpendicular  to 

Pith’  *  <*>  C  “»»•  P*«iug  through  scintillation 

layer  a,  aenith  angle  a,  and  h  .  distance  from  ground  to  scintillation  layer. 


Cloudy  Sky 


‘  xor  scanner. s  with  a  smal’  instantaneous  field  of  view^i 
input  N  (ji  watts)  to  a  detector  viewing  aperture  A  (cm^)  is 


wnere  s  is  scan  rate  (rad/ sec),  fjHfcf^  *s  the 'time  frequency  band-pass  of  the 

system,  W  is  spatial  frequency  (cpr)  and  Q  is  the  function  plotted  above.  This 

evaluation  is  accurate  only  when  the  largest  angular  dimension  of  the  field  is  less 

i 

than  4/s  *  More  complex  geometriei  arc  treated  in  the  literature. 


.mp.  ypi  •#  ’ 


2.  6i  Atmospheric  Scintillation 


Calculations  of  the  scintillation  of  an  extended  source  are 
reported,  in  Appendix  B.  A  summary  of  these  results  will  be  pre- 
senced  here. 

One  may  understand  the  calculated  results  by  means  of  the 

■f  £ 

following  model:  The  refractive  cells  move  as  they  decay  across 
the  projectioh  of  the  telescope  onto  the  turbulent  layer.  Since  most 
scintillation  occurs  near  the  tropopause,  we  take  the  height  of  the 
turbulent  layer  as  the  height  of  the  tropopause,  about  10  km  at 
middle  latitudes.  There  are  various  sizes  of  refract  ve  cells,  and 
which  ones  contribute  the  most  to  scintillation  must  be  determined. 
The  geometry  is  shown  in  Figure  15. 

To  Star 


Telescope 


_ _  Turbulent  Layer 


=  10  km 


Earth 


Figure  15 


Consider  a  turbulent  cell  of  diameter  d.  Its  diffraction  pat¬ 
tern  on  the  ground  has  ait  angular  width  0  of  about  0  —  ^/d,  where 


v, 


A  is  the  wave  length  of  light.  The  angular  extent  of  the  cell  at  the 
telescope,  fi,  is  given  by  fi  -  d/L,  where  L  is  the  height  of  the 
turbulent  layer.  The  geometrical  shadow  of  the  cell  will  be  equal 
to  the  diffraction  pattern  of  th*  cell  when  0  =  />.  For  smaller  cells, 
i.  e.  ,  0  >  )6,  the  geometrical  shadow  essentially  disappears.  For 
larger  cells,  i.  e.  0  <  fb,  the  cell  population  is  much  less  than  for 
the  smaller  cells.  The  diameter  of  the  cell  for  which  0  =  fi  is 
dQ  *  TAL  .  Hence,  cells  of  about  diameter  d  =  dQ,  will  mainly 
affect  scintillation.  As  a  r*5ult,  the  scintillation  of  a  star  as  received 
by  a  point  detector  could  be  considered  as  having  come  through  a  cell 
of  diameter  d  *  dQ  a  "fAI<  •  The  length  dQ/2  is  the  correlation 
length  on  the  ground  for  t\tfo  detectors,  by  a  similar  argument.  An 

X 

extender*  source,  therefore,  may  be  considered  as  beirg  made  up 

of  cells  of  length  d  .  The  scintillation  noise  from  each  cell  is  uncor- 

I 

related.  Those  cells  which  are  not  on  the  edge  of  the  source,  i.e.  , 
internal  cells,  give  no  *..»•  tribution  to  the  scintillation,  due  to  the 
fact  that  light  scatters  one  internal  cell  to  another.  An  extended 

source  projected  onto  the  turbulent  layer  is  shown  in  Figure  16  below. 


Extended  Source 


Typical  Cell 


;ure 


. .  “  uroKen  UP  >™°  cells,  each  of  size  d  ,  twice  the 

correlation  length.  Since  the  scintillation  contribution  from  each 

cell  is  incoherent,  the  total  contribution  is  gfvo'n  by  multiplying  the 
scintillation  noise  by  the  square  root  of  the  number  of  cells.  How¬ 
ever.  the  source  in  each  cell  is  coherent,  so  the  total  source  is 
given  by  multiplying  the  source  in  each  cell  by  the  number  of  cells. 
Hence,  if  the  scintillation  noise-, o-signal  noise  ratio  is  measured 

for  a  star,  one  can  calculate  very  closely,  using  the  equation  given 
below,  the  noise -to- signal  ratir,  _ 


instrument 


where 


«un  a  Noise -to- signal  ratio  of  the 
Number  of  correlation  cells  on  th 
of  the  sun. 


ring  of  correlation  cells 


i  constant  Which  accounts  for  the 
tion  cells’ on  the  sun. 

,  3 

illustrated 


geometry  is 


in  Figure  17  below 
Sun 


Turbulent  Layer 


0  =  l/l  10  radians 


Earth 


-v 


For  purposes  of  the  present  discussion  we  may  con¬ 
sider  the  cells  as  squares  on  the  edge  of  the  sun.  Now  the  num 
ber  of  correlation  cells  on  the  circumference  of  the  sun  is 


becomes 


sun 


r  (2/Vn)(d  /D) 


To  calculate  k,  one  may  observe  that  roughly  half  of  the  scintillation 
noise  from  a  given  cell  will  fall  on  one  of  the  two  adjacent  edge  cells 
and  should  therefore  not  be  counted  as  noise.  This  indicates  that 
k  ~  i/Z.  A  more  exact  calculation  shows  that 


star 


Since  (n/s)fitar  S  Equation  (8)  gives  k  -  i/2.  Now,  the 

cell  size  dQ,  is  equal  to  (i\L)  ^  ,  and  the  diameter  of  the  sun  pro- 

» 

jected  on  the  turbulent  layer  is  D  s  OL.  Substituting  all  these  relations 
into  Equation  (6)  gives 


—  i  .  a,-...,,,,.- _  i s_  i 


E.  H.  PLESSET  ASSOCIATES,  INC, 


3/2^  /t, 3/4 


(n/ s)gun  «  [V(V«)6“a'r~+  1  -  l]  (2/fT)(iiO)  '  (X/L) 


(9) 


The  measured  value  of  (n/s)star  is  approximately  1  for  an  instrument 
whose  aperture  is  less  than  the  area  of  a  correlation  cell  roughly 
Ac  —  (dQ/2)^  ~  13  cm^.  With  this  value  and  with  5000  X, 
Equation  (9)  becomes 

(n/s)  -  1.0  x  10"5 

'  9  sun 

•  i' 

A  more  exact  calculation,  given  in  Appendix  B,  yields  9.  42  x  10  ^ . 

The  stellar  scintillation  decreases  as  the  radius  of  the  telescope, 
for  telescopes  of  larger  diameter  than  the  correlation  length.  This 
fact  may  be  understood  by  noting  that  the  telescope  may  be  broken  up 
into  correlation  areas.  The  contribution  ox  each  area  is  incoherent, 

whereas  the  source  is  coherent.  Hence,  the  total  scintillation  noise- 

? 

to-signal  ratio  is  inversely  proportional  to  the  square  root  of  the  num- 

| 

ber  of  correlation  areas  Within  the  telescope  aperture.  This  is  the 

same  as  to  the  scintillation  noise-to- signal  ratio,  varying  inversely 

jl 

with  the  radius  of  the  telescope. 

'r  » 

2.  62  Frequency  Spectrum  of  Scintillation 

The  frequency  spectrum  of  scintillation  is  dependent  upon  the 
aperture  of  the  telescope  used,  the  size  distribution  of  the  refracting 
cells,  their  lifetime,  theil  height, and  their  velocity  of  motion. 

Figure  18,  a,  b  and  c,  prisents  characteristic  results  of  scintillation 
studies.  For  resolution  of  the  class  of  3  x  10  radians,  seeing  and 
scintillation  are  independent  of  zenith  angle  if  z  <  60°.  For  poor 


••14 *»*#•«  iai|li|tb2 


••l*»|«MW  4»*l»»|<ikJ 


!•»»  '•»!?  (**|*i|o*3 


resolution  (less  then  10'5  radians)  the  variation  is  approximately 

(sec  z)1/2,  and  the  amplitude' of  . scintillation  increases  roughly  in 
the  same  manner. 

The  scintillation  frequency  spectrum  results  from  the  scintil¬ 
lation  of  turbulent  cells  of  different  sizes.  These  cells  move  at  wind 
velocity  and  decay  with  a  lifetime  of  roughly  .  05  sec,  as  was  mentioned 
above.  Hence,  for  any ‘source  of  reasonable  size,  the  turbulent  cells 
move  across  it  before  they  decay.  As  a  result,  each  frequency  in 
the  scintillation  spectrum  is  due  to  a  cell  of  different  size,  with  the 
larger  cells  giving  lower  frequencies. 

The  experimental  data  given  in  Appendix  E  indicates  that  the 
frequency  spectrum  for  ah  extended  source,  in  this  case  the  sun,  is 
roughly  the  same  as  that  for  a  point  source  like  a  star,  except  that 
lower  frequencies  give  relatively  more  scintillation  than  the  higher 
ones  for  the  case  of  an  exlended  source.  This  may  be  understood  in 
the  following  way:  In  Equation  (7)  the  scintillation  noise -to- signal 
ratio  for  an  extended  source  ic  shown  to  be  proportional  to  the  noise- 
to- signal  ratio  for  a  point  source  multiplied  by  the  ratio  of  the  size 
of  a  turbulent  cell  to  the  projected  diameter  of  the  source,  to  the 
3/2  power.  There,  the  cell  used  was  interpreted  as  the  average  cell. 
Now,  since  each  cell  size  contributes  to  a  different  frequency  in  the 
scintillation -power  spectra,  Equation  (7)  may  be  interpreted  as  the 
frequency-dependent  noise-to-signal  ratio  of  the  point  source  and 
extended  source.  It  then  predicts  that  the  scintillation  as  a  function 
of  frequency  for  an  extended  source  has  a  faster  variation  with  fre- 
quency  than  for  a  point  source,  which  seems  to  match  the  data 
obtained  and  reported  in  Appendix  E. 
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i  appendix  a 

Resonan.  Scattering, and  Excitation  Transfer  for  the  Lyman- Birge- 
Hopfield  Series  and  for  the  First  Positive  System  of  Nitrogen 

by  E.  Zipf 

Excitation  of  the  Lyman- Birge -Hopfield  Sytem  of  Nitroeen  by  Resonance 
Scattering  — - - - 1 - - - 

The  Lymah- Birge -Hopfield  system  of  nitrogen  (LBH) 

is  a  conspicuous  feature  in  the  ultraviolet  auroral  spectrum.  The 

system  is  emitted  by  nitrogen  molecules  in  a  radiative  transition 

from  the  a  nj  state  to  this  ground  state  X*2+g  .  Electric  dipole 

transitions  from  the  a1^  state  to  the  ground  state  are  forbidden  by 

the  symmetry  selection  rule  g  -/-»g  so  that  the  a1*  state  is 

g 

meta  stable. 

For  the  most  part  the  brightest  bands  of  the  Lyman- Birge- 
Hopfield  system  are  emitted  in  the  wavelength  region  1100A-2000A; 
strong  absorption  by  O 2  a|  these  wavelengths  precludes  any  possi¬ 
bility  of  observing  these  bands  from  the  ground.  In  auroras  the  LBH 
system  appears  to  be  excited  by  electron  bombardment  as  well  as  by 
cascade  radiation  from  higher  singlet  states.  Presumably  these 
processes  would  also  lead  to  the  emission  of  the  LBH  bands  following 
an  atomic  oomb  explosion.  In  this  connection,  however,  the  following 
question  can  be  raised:  Does  resonance  scattering  of  ultraviolet 
photons  from  the  explosion  by  ground  state  nitrogen  molecules  con¬ 
tribute  in  a  significant  way  to  the  emission  of  the  Lyman-Birge- 
Hopfield  system? 


The  photon  emissiort  rate  per  second  per  N,  molecule  in  the 


v'-v"  band  of  the 
molecules  in  the  v" 


due  to  resonance  scattering  by  N,(XJ2+  ) 

1  g 
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2  A  ,  ,, 

_  tre  .  v',  v" 

8  *  ^avST'o.v'l ^Tv,t 


where  f  ,  is  the  f-value  for  absorption  from  v"  =  0  to  v'  in  the 
o#  v 


upper  state  and  F  p  is  the  incident  flux  of  photons  per  unit  frequency 


interval  per  unit  time  per  unit  area  (Chamberlain,  1961).  The  oscil¬ 


lator  strength  f^(  v„  is  related  to  the  transition  probability  Ay,  v„  by 


f  ,  „  =  1.  5  x  10‘8  X2,  A  ,  „ 

v,v"  v',v"  v,v" 


(1/TV.)  sS  Ayi  vn 


T  ,  is  the  radiative  lifetime  of  the  vibrational  level  v'.  Lichten  (1957) 
v 

-4  ^ 

has  obtained  an  approximate  value  of  1.  7  x  10  sec  for  T(a'Tr  ). 
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Given  the  minimum  signal  detectable  by  a  satellite  LBH- photometer 


one  can  determine  from  Equation  (1)  (making  allowance  for  absorption 


by  0^)  what  class  of  bombs  can  provide  the  required  flux  for  an  observ¬ 
able  signal.  The  prospects  for  a  workable  detection  system  based  on 


the  excitation  of  the  LBH  system  by  resonance  scattering  alone  are  n^t 


very  good  because  the  comparatively  long  lifetime  of  the  a  tt  state 

O 


results  in  a  low  scattering5 efficiency.  A  useful  rule  of  thumb  may  be 
derived  from  the  recent  dayglow  experiments  of  Fastie  et  al.  (1964), 
who  found  that  under  normal  conditions  (no  auroral  activity)  there  is 


no  observable  LBH  emission  (<100  rayleighs)  in  the  dayglow.  This 


would  imply  that  the  flux  incident  on  the  upper  atmosphere  in  the  11Q0A- 


2000A  region  of  th«  spectrum  from  a  bomb  explosion  would  have  to  be 


considerably  more  intense  than  the  solar  flux  at  these  same  wave 


lengths  to  produce  a  detectable  effect. 
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II.  Excitation,  of  the  First  Positive  System  of  Nitrogen  oy  Collisions 
of  the  Second  Kind  with  Metastal  le  Nitrogen  Molecules  in  the 


a  it  State 

- g_ — — 


Oldenberg  (1959)  has  advanced  the  hypothesis  that  the  reaction 
N2(a%  )  +  N2 -^N2(B\  )  +  N 


(4) 


is  the  principal  excitation  mechanism  for  the  First  Positive  bands  of 

.  ...  ' 

nitrogen  in  the  Thompson- Williams  (1934)  afterglow.  This  afterglow 

O 

is  produced  by  bombarding  N2  at  a  pressure  of  3  x  10  mm  Hg  with 
electrons  at  energies  greater  than  10  ev.  The  experimental  data  are  in 
qualitative  agreement  with  Oldenberg1  s  theory.  However,  no  esti- 

■  it  •  i,< 

are  available  of  the  magnitude  of  the  cross  section  for  process 

rJ  '  rj'k 

(4).  More  experimental: work  is  desirable. 
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The  Effects  of  Overcasts  on  the  Light  Received  by  the  VELA  System 

by  J.  DeGroot 
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This  appendix  has  been  abstracted  from  a  series  of  reports 
prepared  at  E.  H.  Plesset  Associates  under  AFTAC  Contract 
No.  AF  33(657)- 13599.  These  reports  discuss  the  effects  of  clouds 
and  overcasts  on  the  signature  of  X-ray  induced  fluorescent  light 
received  by  the  VELA  system.  The  source  of  the  X-rays  is  a 
nuclear  burst  in  space.  Specifically,  the  results  of  this  work  are 
presented  in  the  form  of  curves  of  the  intensity  of  light  received  by 
a  VELA  detector  with  an  overcast  present  relative  to  the  intensity 
yecieved  without  an  overcast.  The  nuclear  burst  is  assumed  to  occur 

at  the  zenith  of  the  detector  and  far  enough  away  so  that  the  X-ray 

| 

pulse  is  a  plane  wave  when  it  strikes  the  fluorescent  layer.  The 
fluorescent  layer  is  arsumed  to  be  optically  thin,  and  the  fluorescence 
is  assumed  to  occur  in  a  very  thin  layer  at  an  altitude  of  80  km. 

With  these  assumptions  for  the  source  of  light,  a  uniform  over- 

I  V 

cast.infinite  in  two  directions,  was  used  in  the  calculations.  The 


geometry  is  shown  in  Figdre  1  below: 
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The'time  dependent  transmission  of  the  overcast  was  calculated 
using  diffusion  theory.  The  detector  was  assumed  to  have  an  iso 
tropic  response  out  to  80°  from  the  vertical. 


The  calculat?ons  were  done  on  an  IBM  1620  binary  digital 
computer.  Overcast  thickness,  height,  and  diffuse  transmission 


were  parameters.  Calculations  were  made  with  overcast  thick 


2  km,  4  km,  and  6  km,  and  with  overcast  heights 
4  km.  and  6  km.  Overcast  diffuse  transmissions 


which  were  used,  are:  .07,  .15,  .  25,  .  35,  and  .  50 


The  resulting  set  of  curves  are  shown  in  Figures  2-11,  For 
each  curve,  the  relative  intensity  received  by  the  VELA  system  is 
shown  without  clouds  as  a  f\  notion  of  time.  The  origin  in  time  is 

when  the  first  light  is  received  at  the  detector.  The  curves  are 

-  } 

normalized  to  the  light  intensity  received  without  an  overcast  at 

n  , 

the  origin  in  time.  The  Overcast  thickness  and  height  are  constant 
on  each  figure,  while  the  overcast  diifufe  transmission  is  varied. 
One  effect  that  can  be  seen  in  all  the  figures  is  that  the  over- 

.7  j  ? 

cast  can  sharpen  the  light  pulse.  This  occurs  for  a  thin  overcast, 
i.  e. ,  T^qO  —  .  5  .  Ths  physical  reason  for  this  is  that  the  over- 
ast  acts  only  to  redirect,  not  to  diffuse,  light.  Hence,  the  over- 

cast  may  be  considered  to^be  a  i\ a,  detector,  which  has  an  extra 

I 

cosine  factor  in  the  light  intensity  received  over  that  received 

» 

£ 

without  clouds. 
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On  Scintillation  and  Atmospheric  Turbulence 
by  G.  Milne 

'■  Summation  o£  apparently  small  light  sources  (stars  or  distant 
terrestrial  sources)  results  from  optical  path  variations  apparently 
due  primarily  to  thermal  inhomogeneity  of  the  intervening  medium. 

For  horizontal  paths  these  inhomogeneities  are  normally  distributed 
more  or  less  uniformly  along  the  line  of  sight.and  the  resultant  degrad, 
ation  of  optical  imagery  may  be  dealt  with  fairly  well  on  the  basis  of 
statistical  models.  For  vertical  paths  the  situation  is  rather  different, 
the  most  significant  difference  undoubtedly  being  the  considerable 
degree  of  stratification  which  is  normally  present.  The  various  air 
layers  in  the  line  of  sight  are  not  reaaily  accessible,  so  that  their 
properties  must  be  inferred  from  observations  which  must  usually  be 
made  on  the  ground.  Assuming  that  at  least  a  large  part  of  the  scintil¬ 
lation  observed  in  the  case  of  a  star  arises  from  a  single  disturbed 
layer  at  a  considerable  altitude,  the  velocity  (including  speed  and 
direction)  of  such  a  layer  may  in  principle  be  inferred  from  correla¬ 
tion  of  observations  with  two  telescopes  separated  by  a  Known  distance. 
Further,  the  apparent  angular  velocity  of  such  a  layer  may  be  deter¬ 
mined  from  correlation  of  observations  of  two  stars  of  known  angular 
separation,  from  which,  knowing  the  linear  velocity,  the  altitude  of 
the  layer  may  be  found.  In  practice  these  observations  turn  out  to  be 
very  difficult,  due  in  part  to  signal-to-noise  problems  associated  with 
the  low  radiant  energy  density  available  from  stars,  but  primarily 
due  to  the  nature  of  this  tui-bulence,  since  the  lifetimes  of  individual 

segments  of  the  spatial  turbulence  pattern  are  of  the  same  order  as 

-56- 


the  time  required  for  them  ,o  move  their  diameter,  so  that  quite 
sophisticated  correlation  techniques  are  required,  with  long  times 
°f  Observation.  Reasonable  progress  has  nevertheless  been  made 
in  this  direction.  To  my  knowledge,  no  actual  measurements  of 
apparent  angular  velocity  have  been  made.  I,  has  long  been  observed 
that  double  stars  scintillate  together,  indicating  that  for  separations 
of  the  order  of  a  few  seconds  of  arc  any  correlation  time  interval 
must  be  very  small,  whieh  is  in  agreement  with  other  crude  observa¬ 
tions  of  my  own,  mentioned  below.  The  distribution  of  adequately 

bright  stars  is  not  favorable  for  such  observations,  so  that  on  the  whole 

the  usr  of  stars  for  atmospheric  tu’*bulp>nrA 

^  c  [U‘  ouience  measurements  seems  of 

rather  limited  value. 

2-  There  still  seems  to  be  a  need  for  a  great  deal  more  observation 
of  scintillation  for  vertical  paths,  with  an  effort  to  relate  these  to  cur¬ 
rent  theories  on  turbulence,  and  perhaps  relate  1,  quantitatively  with 
variations  in  temperature  or  other  quantities.  Probably  the  most 
significant  contribution  that  could  be  made  as  an  initial  step  would  be 
the  determination  of  velocity  and  altitude  for  turbulent  layers  and 
possible  association  of  these  with  regions  of  shear  such  as  the  tropo- 
pause.  In  this  connection,  further  consideration  might  be  given  to 
variations  of  longer  period,  of  the  order  of  0.  1  second  to  several 
seconds,  and  one  means  of  approaching  this  might  be  by  direct  observa¬ 
tion  of  schlieren  effects  on  the  limb  of  the  sun,  moon  or  planets.  While 
such  observations  are  strictly  not  identical  with  scintillation  observa¬ 
tions,  they  are  similar,  and  would  appear  to  provide  a  valid  approach 


4ft  * 


to  the  obtaining  of  data  on  physical  conditions  of  the  atmosphere 
overhead.  Rough  estimates  of  the  velocity  of  atmospheric  disturbances 
along  the  limb  of  the  sun  and  moon  indicate  values  of  angular  velocity 
in  the  order  of  several  minutes  of  arc  per  second,  and  in  the  case  of 
the  moon  some  correlation  has  been  observed  between  this  angular 
velocity  and  movements  of  shadow  bands  for  stars  in  the  same  part 
of  the  sky.  Observations  might  be  made  either  photo-electrically 
or  photographically,  and,  with  appropriate  correlations,  the  opportunity 
exists  for  obtaining  both  velocity  and  altitude  f  turbulent  layers. 
Furthermore,  in  addition  to  the  advantage  of  greater  available  energy, 
there  is  also  the  possibility  of  obtaining  data  during  the  daytime  as 
well  as  at  night,  which  would  be  very  difficult  if  only  stars  are  used. 
Since  the  time  resolution  required  may  severely  limit  exposure  time 
m  photographic  observations,  it  is  possible  the  use  of  image  orthicons 
might  be  valuable  here.  Such  observations  would  include  simultaneous 
viewing  with  two  separated  apertures. 

3.  There  may  be  interesting  possibilities  in  correlating  results  of 
vertical  ranging  with  observations  of  celestial  objects  near  the  horizon 
since  evidence  of  atmospheric  stratification  may  be  frequently  seen 
here.  Such  observations  might  be  made  on  the  sun.  moon,  planets, 
and  even  satellites  such  as  Echo  I  and  II,  but  knowledge  of  atmospheric 
refraction  may  not  be  adequate  for  obtaining  quantitative  information. 

4.  While  interest  in  stellar  scintillation  has  come  about  chiefly  in 
connection  with  astronomical  observation,  and  more  recently  m  con¬ 
nection  with  "seeing"  Whigh  altitude,  the  most  important 
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significance  of  turbulence  observations  might  well  lie  ;n  the  field 
of  meteorology,  particularly  if  reasonably  simplified  observation 
techniques  can  be  developed  which  would  permit  simultaneous  observa¬ 
tions  from  a  number  of  stations.  It  should  be  possible  to  follow  such 
Phenomena  as  jet  streams,  and  perhaps  more  subtle  fluctuations  m  air 
mass  movements.  The  limitations  of  optical  wave  lengths  under  cloud 
cover  would  make  it  desirable  to  work  at  other  "windows"  at  longer 
wave  lengths,  but  this  would  undoubtedly  be  very  difficult. 

5.  In  connection  with  (2),  it  might.be  of  considerable  interest  to 
cheok  correlations  with  time-lapse  photography  (or  similar  observa¬ 
tions)  on  light  cloud  cover  such  as  cirrus,  altocumulus,  and  perhaps 

even  jet  contrails.  Strong  turbulence  may  frequently  be  seen  associated 
With  scattered  cumulus-type  cloud  cover  at  night. 


E.  H  PLE5SET  ASSOCIATES,  INC. 


r'N 

^r\ 

Ct 


APPENDIX  D 

ATMOSPHERIC  SCINTILLATION  OF  AN  EXTENDED  SOURCE 
by  Wendell  A.  Horning  and  J.  J.  McClure 
1.  Statement  of  Problem 


A  detector  is  considered  to  view  the  sun  or  a  planet 
through  the  atmosphere.  The  detector  output  is  taken  to  be  propor¬ 
tional  to  the  integral  of  light  intensity  over  the  detector  area.  Even 
if  the  light  source  is  steady,  detector  output  will  vary  randomly  as  a 
result  of  atmospheric  turbulence.  It  is  desired  to  calculate  the  vari¬ 
ance  Ox  the  fluctuations  in  this  output,  their  spectral  density,  as  well 
as  their  correlation  length.  If  two  pencils  of  light  cross  the  turbulent 
layer  at  two  positions  separated  by  a  distance  q  smaller  than  the  scale 
length  w  of  the  turbulence,  they  suffer  amplitude  and  intensity  changes 
which  are  correlated,  while  there  may  be  negligible  correlation  if 
q  »  w  .  Thus,,  there  is  a  statistical  correlation  between  the  outputs 
of  two  detectors,  to  be  calculated  as  a  function  of  their  separation,  b. 

It  will  turn  out  that  fluctuations  in  detector  output  are  sensitive 
to  both  the  size  and  the  shape  of  an  extended  source.  A  description  will 
be  given  of  a  method  for  deducing  at  least  the  gross  features  of  source 
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size  and  shape  from  an  analysis  of  these  fluctuatxons.  A  second 
use  for  knowledge  of  turbulence  produced  fluctuations  is  in  the  design 
of  experiments  to  measure  the  altitude  distribution  of  turbulence, 
which  is  at  present  not  known  in  detail.  More  knowledge  of  this 
distribution  would  contribute  to  the  safety  of  jet  aircraft  and  would 

be  of  use  in  studies  of  weather,  as  well  as  in  the  design  of  optical 
systems. 


2'  Two  Poi^t  Sources  and  Two  Point  Detectors 

The  problem  of  Section  1  can  be  approached  through  a  simpler 
problem,  as  illustrated  in  Figure  Z.  Sj  and  S.,  are  two  point  sources 
of  light  separated  by  distance  a,  and  at  distance  zg  from  the  turbu¬ 
lent  layer.  and  D2  are  point  detectors  separated  by  distance  b. 
and  at  a  distance  zr  from  this  layer.  Light  from  S.  may  reach  D^by 
the  path  S1D1  to  produce  a  displacement  in  the  optical  field  of  the  form 
AilC  11  ,  where  and  are  an  amplitude  and  a  phase.  The 
total  displacement  at  resulting  from  the  two  sources  will  be 

A11eiSU-+  A12eiS12 
The  similar  displacement  at  D2  will  be 


(1) 


$2  c  A21®lS2i  + 


A22*iS22 


(2) 


A  path  such  as  Si may  be  formally  called  a  ray,  but  we  do  not 
rasan  to  exclude  effects  of  diffraction. 

The  power  reached  at  will  be  ,  where  a  star  indi¬ 

cates  a  complex  conjugate.  Or, 

i 

J1  *  AU  *  A122+  2AllAUcos<5irSi2>  (3) 

If  tne  light  sources  Si  and  S2  are  incoherent,  as  is  assumed,  then  the 

phases  and  S12  will  be  Jmcorrelated,  and  the  cosine  term  in  «) 

may  then  be  ignored.  Intensities  at  a  single  detector  will  now  be  addi¬ 
tive,  or 
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J1  =  J11  +  J12  .  "2  =  J21  +  J22  , 

where  J  for  example,  is  an  abbreviation  for  A,,*.  Correlations 

between  two  quantities  such  as  Ju  or  or  between  J, ,  and  J 

may  not,  however,  be  ignored.  The  essence  of  the  nrob’  21- 

u  01  tnc  problem  consists 

kln2  S1ICh  '«~»«‘on.  into  account  in  the  proper  way. 

Let  j  and  j.,  be  the  fluctuations  in  the  output  of  detectors  1  and 

o-od^e  ‘b  r  T  flUCtU2.ti0n  “  ttat  «*•«*“  of  the  output  of  detector  l 
P-  uce  by  Ught  from  source  1.  as  shown  in  Figure  2,  and  so  on. 
x.nen, 

ji  e  J'n  +  hz 


.2  .2  ,2 

Ji  =  J11  +  J 12  + 


Write 


2  _ 

'1  =  J1 


_  2  _  — 
'  fll  =Jil 


2  _  “  2 
*12  =Ji2 


Lll,  12  ~  «?‘l  1J 12 


Ih  (6)  the  bar  indicates  a  statistical  average  which  may  be  taken  as  a 
ime  average  since  the  random  properties  of  the  turbulence  are 
assumed  to  be  stochastic.  The  symbol  a  indicates  the  variance  of 
some  quantity,  a,,  for  example,  is  the  variance  or  runs  value  of 
■he  fluctuations  in  the  output  of  detector  1.  ,  -j.  thc  variancc  „f 

‘  P°rtl0n  *e  °Mput  °£  *>“«<>*  1  which  arises  from  source  1 
and  ,0  on.  The  symbol  E  in  (6)  indicates  thc  average  of  the  oroduct 

of  two  fluctuations,  and  is  called  the  correlation  between  these  flue- 
tuations.  By  (5)  and  (6) 

"T2  *  “'ll  +  °rf^2+  2RU,12  (7 

which  is  the  rule  for  finding  the  variance  c  ,  of  the  sum  of  the  two 
uctuations  J,  j  and  j12  which  may  be  correlated.  The  correlation 


l; 


•b.  ■  f 

„  . 
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between  the  total  fluctuation,  of  actors  1  and  2,  in  the  above 
notation,  will  be 


© 


© 


12 


h>z  *  zf  ■'  'Km  , 


V 


(8) 


a  summation  which  is  easily  exto-ied  to  the  case  of  a  source  composed 
o  many  discrete  points.  It  is  customary  to  define  the  correlation 
coefficient/)  between  two  quantities,  say  j  n  and  ■"  y  by  the  equation 


At. 


12  *  Ru,  tz^  u"  iz 


(9) 


This  dimensionless  quantity/)  must  lie  between  -  1  and  +  1.  If  ,he 

•arious  point  sources  are  of  equal  strength  and  if  the  turbulence  is 
statistically  uniform  in  space,  then  o  =  o  .  ,  and  (9) 

may  be  written  14 


r12A2 


■sy* 

/Ay 


\M,2V 


(10) 


Extended  Source 


If  the  source  covers  4n  area  A,  but  the  detector  at  D.  has  a 
diameter  small  compared  to  the  scale  length  of  the  turbulence,  then  the 
to.al  fluctuation  of  the  output  of  D,  is  given  by  the  new  formula 


j]  -  ds 


UD 

‘  indicatls  ?fsition  on  «“  «•»«.,  and  d?  is  an  element  of  source  area  J 
,  Jl(l)  dS  1S  *  diStinct  random  ™*bU  for  each  distinct  element 
of  area  ds.  These  variables  are  infinite  in  number  and  correlated. 

Their  joint  probability  distribution  is  hence  infinitely  dimensional, 
ut  the  usual  statistical  averages  may  still  be  taken.  The  square  of 

jl  may  be  expressed  by  the  double  integral 

£ 


M  Z 
-1  " 


* 

i(s)j ^s^ds  ds' 


(12) 
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The  variance  of  j1  will  then  be 

ffl  ^1  =  '/a.J A  Rl^s* 


s')  ds  ds' 


(13) 


where  R^s,  s')  ds  ds'is  the  unnormalized  correlation  between  those 
two  contributions  to  fluctuations  at  the  detector  Dj  which  arise  from 
area  elements  ds  and  ds'.  The  form  of  the  function  R^s,  s')  is 
deducible  from  Tatarski.  If  the  source  area  is  of  uniform  brightness 
and  if  the  source  diameter  subtends  only  a  small  fraction  of  ^  adian 
at  the  det  ctor,  then  is  a  function  of  |s-s'J  .  If  the  limb  darkening 
of  the  sun  is  to  be  treated,  then  R1  may  be  written  approximately 

*C(s)Ru(|s-J'|)  (i4j 

where  R^  is  the  form  for  Rj  when  the  limb  darkening  is  ignored  and 
G(£),  unity  at  the’ sun's  center,  is  the  fractional  darkening  as  the 
source  point  moves  toward  the  sun's  edge.  Equation  (14)  is  quite  good 
for  the  sun,  because  in  the  symbols  of  Figure  2, 

Ds/zo  <<  w/zr 

where  Dg  is  the  sun's  diameter  and  w  =  10  cm  is  the  scale  length  for 
atmospheric  turbulence.  That  i-.,  within  a  correlation  length  pro¬ 
jected  *o  the  sun,  its  brightness  changes  by  a  very  small  fraction. 

Tc  compute  the  correlation  R12  between  the  outputs  of  the  two 
poirt  detectors,  we  write  for  these  outputs 


j-  =o/ji(s)ds,  j2  ,fhCs)  ds 


(15) 


For  each  infinitesimal  pencil  of  light  from  source  to  detector  there 
is  a  distinct  random  variable.  Two  pencils  a  he  same  only  if  they 
start  from  the  seme  source  area  element  and  end  cn  the  same  detecto; 
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The  result  is 


R12  =  R12^s'  ®') 


in  whir  R12{s,8')  ds  ds'  is  the  unnormalized  correlation  between 
increments  of  the  detector  power  which  may  or  may  not  arise  from 
the  same  source  area  element,  but  which  must  be  delivered  to  dis¬ 
tinct  detectors.  As  in  the  similar  result  (13),  the  function  R 
depends  only  on  s-s'  if  the  source  is  uniform  in  brightness.  If  not, 
a  limb  darkening  function  may  usually  be  introduced  as  in  (14). 

If  the  two  detectors  at  Dj  and  D£  ,  as  well  as  the  source,  extend 
over  significant  areas, ‘then  the  fluctuation  of  the  output  from  an  area 
element  dr{  of  Dj  is  j^Fj)  dFj  .  The  total  fluctuation  at  Dj  will  be 

d*  *  where  Bj  is  the  sensitive  area  of  Dj.  The  variance 
of  jj  will  be  <r  ^  ,  given  by 

°iz  *  Ads-4d6'-4Idi'/ B1dr' (n 

where  R(r,  s,  r',  s')  dr  ds  dr'  ds'  is  the  unnormalized  correlation 
between  an  element  of  the  power  to  Dj  which  flows  from  dF  to  dF  and 
a  second  such  element  which  arises  from  dF'  and  is  delivered  to  dF'. 

*2  *  WhlCh  &pplies  ^  the  detector  D2,  is  given  by  (17)  after  the  sub¬ 
script  1  is  replaced  by  2.  If  the  source  is  of  uniform  brightness  over 
its  area  and  if  both  source  and  detector  diameters  subtend  only  a 
small  traction,  of  a  radian  at  the  turbulent  layer,  then  the  function  R 
of  (17)  can  be  simplified  to  the  form 


R  ~  R(f  r-r'  | ,  j  s-s'j  ,  f) 

wnere  $  is  the  azimuth  angle  between  the  vectors  F-F'  and  s-s'  .  Limb 
darkening  may  be  introduced  as  before. 

A  derivation  similar  to  that  for  (17)  shows  that  the  correlation 
between  outputs  of  the  two  extended  detectors  D,  and  D  "is*'3 
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©I 


*  ^  12  “  dr/  dp  R(r,s.P,s') 


where  Bj  and  B2  are  the  sensitive  areas  of  the  two  detectors. 


4*  Correlation  Between  Two  Narrow  Beams 


The  central  diffraction  maximum  produced  by  an  eddy  of  size  w 
will  have  an  angular  half  width  oU/w.  where  ^  is  the  wave  length  of 
light,  this  from  elementary  principles.  Suppose  w  is  so  small  that 
A/w  is  large  compared  to  the  angle  subtended  by  w  at  the  detector 


on  the  ground. Then,  where  zy  is  altitude  of  the  turbulence. 


A/w  »  w/z 


If  (19)  is  satisfied  by  w,  then  the  diffraction  patterns  of  most  eddies 
will  overlap  at  the  point  detector,  and  there  will  be  statistical  can¬ 
cellation  of  their  fluctuations  at  the  detector  even  though  this  cancel¬ 
lation  will  not  be  complete.  If,  on  the  other  hand. 


A/-W  <<  w/z 


then  at  any  given  time  the  diffraction  pattern  of  only  one  eddy  of  size 
w  will  significantly  affect  the  detector,  and  small  changes  in  its  cen¬ 
tral  position  will  but  little  affect  the  fluctuations  of  detector  output. 
Kence,  intuitively,  eddy  sizes  which  satisfy 


A/w  -  w/z  ,  w  ~  (/z  )' 
1  r 


will  make  the  chief  contribution  to  detector  fluctuations,  a  result  which 
agrees  with  a  more  rigorous  discussion.  This  assumes,  of  course, 

-t  zr  is  Sufficiently  large  so  that  .he  detector  is  in  the  diffraction 
zone,  not  the  geometrical  optics  zone.  In  practice,  this  means,  that 
zr  is  greater  than  a  few  hundred  meters,  which  includes  the  cases  of 

present  interest.  Ker.ce,  the  correlation  length  of  fluctuations  .at  the  -round 
will  be  of  the  order  of  z*)1'2,  . 


m 


ns 


■  v  1 
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C01reUti0n  “"OlcUntfM  for  fluctuations  at  the 
-  n  an  .tom  a  point  source  is  given  by  Tatarski  on  p.  154  (Ref.  1)  for 

to' *1  n't Tr  ^  tUrbUle“e  ^  °£  Uai£°rm  “tensity  *»«  the  ground 
n  altitude  L.  It  may  her  fitted  quite  well  by  a  function  of  form 

A(r)  rr  r-  *.2/.  2.  9.9 


pM  =  Cjexpl-rVuj2)  .  c2exp(-r2/u22) 


C1  >  c2  •  ui  <  u2 


It  has  the  important  property 


fr 


(r)  dr  =  0 


where  the  integration  is  over  ali  the  horizontal  plane  through  the 
detector.  Equation  (22)  is  equivalent  to  the  assumption  that  the 
turb  t  atmosphere  is  non-absorbing  and  is  a  sufficiently  weak 

fan  s'"  “  ^  aU  ^  WWCh  £aUS  °"  “  £r°m  ab°-  ultimately 
falls  somewhere  on  the  ground.  The  effect  of  the  turbulence  is  then 

°oly  to  redistribute  the  light  on  the  ground.  With  this  assumption 


dr  s  0 


TZr  Z  th'  fluctuation  as  a  function  of  ground  position 

of  T  T  'aV  6°Ur- ' ’■  ^  de£inili°n  16  tha  average 

has  statistic ^  ’  W“er'  8r°Un<i  p°sition-  If  the  turbulence 

theo  T  Pr0P"'I6S  Unchan»i^  ^  time,  then  by  the  ergodic 

theorem  we  have  ° 

/°(r)  “  ^(r1)j(ri+r)  d?4 

.  «  * 

where  the  integration  is  again  over  the  ground  plane.  Integrating 

o (22>-  * £oMows  "°m  w  ~ 

# 1 . 

Clni^  “  C9U9^  a  0 

1  1  2  2  t  /  9 1 


•j 

A  ^  ^ 
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/ 
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A  form  of  the  intensity  correlation  coefficient  more  general 
than  (21)  is 


/)(q)  =  Cjexpf-q2/^2)  -  exp(-q2/u22)  ( 

where  /He,  is  the  correlation  between  the  intensity  from  a  source 
Si  at  a  detector  Dj  and  the  intensity  from  a  source  S2  at  a  detector 
D2*  and  Where  is  the  distance  between  the  two  intersections  of 
lines  SjOj  and  S2D2  with  a  thin  layer  of  turbulence  (see  Figure  2). 
Tatar  sk*  s  plot  of  p  applies  strictly  to  a  uniform  layer  of  turbulence 
extending  from  the  ground  to  a  height  L.  For  the  present  we  con- 
sider  it  to  hold  also  for  a  thin  layer  at  height  L/2.  Then  approximate 
values  of  the  constants  of  (26)  may  be  obtained  by  fitting  his  curve. 

A  more  complete  derivation  of  the  constants  of  (26)  along  with  a 
fit  to  p  possibly  containing  more  than  two  terms  will  be  given  later. 

i 

5*  j^arce  Extended  in  One  Din-mnsir^ 


As  a  helpful  preliminary  problem,  consider  a  line  source  of 
light  2x0  at  distance  from  a  detector  D,as  in  Figure  3.  Light 
from  the  source  traverses  a  thin  layer  of  turbulence  at  height  z  . 
According  to  Sections  3  and  4,  the  variance  <r  in  the  detector  output 
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will  be  given  by 

,2  =  /x°a*/ 

~/-x  v/-x 


|t:i  «xp[-(x-x')2/w12] 

-  c2  exp  |-(x-x')2/w2;2  J 
1  S.  ui  zo/zr  '  w2  “  u2  zo/zr 


i 

) 


;re 


Since  the  turbulence  only  redistributes  light  on  the  ground,  the: 
will  be  a  restriction  on  the  correlation  coefficient  similar  to  (25), 
which  in  the  present  two-dimensional  problem  becomes 


(27) 


r>s 


c\ 

J 


ciul  "  C2U2  =  0 


(28) 


The  first  integration  of  the  first  term  of  the  integrand  of  (27) 


yields 


o  r 


x  -x1  x  +  x*  7 

erf(  - )  +  erf(  — — - ) 


w 


1 


w 


1 


dx1 


(29) 


where 


/“a  _  2 

erf  ( z )  *  £  e"x  dx  ,  z  <  Q 
After  the  integration  over  x1  is  performed,  expression  (29)  becomes 


2x 

+  4x  erf  (  — -  )  l 

o  wi  J 


(30) 


The  second  term  of  the  integrand  of  (27)  yields  a  similar  expression. 

If  the  length  2x0  subtends  a  half  degree  at  the  ground,  as  does 
the  diameter  of  the  sun,  and  if  w^  and  w^  have  sizes  appropriate  to 
actual  air  turbulence,  theh  2x :/wi  and  2x  /w,  are  both  of  the  order 

of  10  .  Then  the  exponential  terms  in  or4  are  negligible.  Then  the  two 

:  fill  :t  i 

* 

* 
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erf  terms  are  nearly  equal  in  absolute  magnitude  and  cancel  with 
the  aid  of  (28).  The  result  is 


The  rms  fluctuation  or  at  a  point  detector  from  a  sufficiently  extended 
source  is  hence  independent  of  source  size  xq.  The  basic  reason  is 
that  the  turbulence  only  redistributes  light  from  the  source.  The 
fluctuations  hence  result  from  the  source  ends  and  represent  an  end 
effect.  We  will  see  later- that  fluctuations  from  an  extended  two 
dimensional  source,  similarly, represent  an  edge  effect. 
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6*  Source  Extended  in  Two  Dimensions 


Now  consider  the  line  source  of  Figure  3  replaced  by  a  rectangle 

^  S‘dCB  2xo 2V0  ■  ^  a  Plane  parallel  ,o  the  layc^r 

o  tur  ulence.  If  zq  »  „r  ,  then  the  problem  will  be  essentially 
unchanged  if  the  source  lies  on  any  curved  surface  area,  provided  the 
normal  projection  of  the  area  on  a  plane  parallel  to  the  turbulence  is  a 

rectangle .with  sides  2*q  and  2yo.  The  variance  <r  at  a  point  detector 
will  now  be  given  by 

<r2  “  Jdxdx'  dy  dy'  jcj  exp[-[(x-x')2  +  (y-y'J^/wj2] 

(32) 

-  c2  exp[  "  (<*-*'>2  +  (y-y')^  w22J 

Where  the  range  of  integration  covers  the  source  area  twice.  Because 

of  the  rectangular  source  shape,  each  of  the  two  terms  of  (32)  factors 

...  >  a  product  of  two  expressions,  each  of  the  form  of  (3C).  Hence 
(32)  becomes  ' 


Y 

"  *  C1W1?  ,2irl/  \  -  Wj)  (2v‘/2yo  .  W[) 

-  c2w22  (2v1/2xo  .  w2)  (2v‘/2yo  .  w  , 


<3i 


which  again  assume,  x0  and  yo  to  be  very  large  compared  to  w,  and  w 
so  that  the  source  dimensions  may  simulate  those  of  an  object  such  ‘ 
as  the  sun.  Because  of  (25),  the  terms  of  (33)  proportional  to  x  y  wil 

cancel,  and  the  terms  linear  in  x  and  v  ^,-n  v*  j  •  .  °  ° 

become8  ear  in  xq  and  yQ  will  be  dominant.  (33)  hence 


<r2  »  2ir1/2(xo+yo)(c2w23  -  c^3) 


(34) 


which  is  proportional  to  source  perimeter  and  hence  exhibits  the 
variance  as  an  effect  characteristic  of  the  source  edge. 


1  4 
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The  result  (34)  can  be  put  into  a  form  consistent  with  intuition 
by  writing 


3  3.  3 

C2W2  ~ClWi  ~cw 


(35) 


where  c  is  new  constant  of  the  order  of  c,  and  w  is  of  the  order  of 

r  " 

the  scale  of  turbulence  increased  by  the  ratio  z  /z  .  w  ic  hence  the 

o  r 

projection  from  the  turbulent  layer  to  the  sun  of  the  scale  of  turbu¬ 
lence,  the  projection  taken  with  respect  to  the  point  detector. 

N  «  2(s0+y0)/w  is  the  greatest  number  of  non- overlapping  areas  A 
of  size  w^  which  can  simultaneously  touch  the  edge  of  the  source. 
Hence,  (34)  can  be  put 


>3 


<r2  =  NA2c 


or  Ac 


1/2 


(36) 


.1/2 


The  quantity  Ae4/  b  is  the  variance  from  a  portion  of  the  source  of 
area  small  enough  6o  that  its  various  parts  are  well  correlated.  The 
factor VN~  then  arises  naturally  from  the  usual  rule  for  adding  the 
effect  of  statistically  independent  fluctuations.  Fluctuations  at  the 
detector  hence  arise  principally  from  a  strip  of  width  w  which  borders 
the  perimeter  of  the  source,  provided  the  source  dimensions  are 
large  compared  to  w. 


£ 

For  the  case  of  a  round  source,  say  the  sun,  (34)  is  replaced 


by 


<r 2  =  2tr3/2Rj[c2W£4^'qWj3) 


(37) 


where  Rg  is  the  source  ra|lius. 

I 


■I 
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^ ’  Correlation  with  an  Extended  Source 

The  aim  of  this  section  is  to  compute  the  correlation  cr 

between  fluctuations  of  the  outputs  of  two  point  detectors  on  the 

g  round , s  e  pa  rate  d  by  a  distance  b.  The  source,  as  in  the  preceding 

section,  is  a  rectangle  with  sides  2xq  and  2yQ.  For  the  present  the 

line  o*  centers  of  the  detectors  is  taken  parallel  to  the  side  of  the 

rectangle  of  length  2x  .  Then 

o 

12  sjdxdydx'  dy^  c^exp  -{(x-x'+x)2  +  (y-y')2}/*^' 

c2  exp[“{(x-x'+x)2.+  (y-y')2}/w„ 

is  the  desired  quantity,  where  the  meaning  of  lengths  and  Zg  is 

given  in  Figure  2.  The  integral  over  x  of  the  first  term  of  the^ntegrand 
of  (38)  yields  ? 


(38) 


cw, 


Vx'+*  -x-x'+x  | 

erf(— ...  j,-)  +  erf  (—£-■.  ■■  )j 


(39) 


The  integral  of  (39)  with  respect  to  x'  and  over  the  range  -x  <  x'  <  x 
yields  °  ° 


Jix  e  W1  f’wl  exp[-(x/wi)2J  -  2x  er^x/wj) 


+  (2xo+x)(tt1/2/2)+  (2xo-x)e>5>r^--£—  j2]l(40) 


K 


+  (w,/2)  exp 


2x  -x 


L 


w 


1 


Similarly,  the  integral  of  the  first  term  of  the  integrand  over  y  and 
y',  both  over  the  range  -yQ  ib  +yo#  yields 

I 

*iy  *  Wj  (-Wl  4  2Wl/2yo) 


J 


(41) 
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which  is  what  (40)  becomes  whoa  x  is  put  equal  to  zero  and  x  is 
replaced  by  v  Hence  (33),  f»Uy  evaluated,  is 

”"12  “  CiIlx^ly  "  C2^2x^2y 
In  (42)  the  te  rms  proportional  to  x  y 

-  .  °  °  aH  cancel  because  of  (25).  If 

x  is  very  small  compared  to  xo  and  yo  .  then  the  terms  linear  in  x 

and  yo  dominate.  But  in  jeneral, terms  linear  in  the  three  quan,itie°s 

*?'  y°.  “ ,  *  arC  0(  U’e  6ime  ordar  cven  ‘°r  *  very  large  source.  A 
f^the^cha rac ter  shown  in  Figure  4 


Figure  4 

•  •  • 

If  the  line  of  centers  of  the  two  detectors,  contrary  to  what 
wae  assumed  above,  make,  any  significant  angle  with  the  sides  of 
the  rectangular  source,  then  the  two  integrals  I  and  I.  become 
Of  symmetrical  form  and  the  plot  of  then  appear,  a  eYn  Figure  5. 


■>  \  ■»: 


Figure  5 


Hence,  the  character  of  <r  12  as  a  function  of  x  is  sensitive  to  the 
angle  between  the  line  of  centers  of  the  detectors  and  the  sides  of 
the  rectangular  source.  This  can  be  interpreted  by  considering  the 
projection  of  the  source  on  the  layer  of  turbulence,  first  with  respect 
to  detector  and  with  respect  to  D^.  If  these  projections  are  as 
in  Figure  6  a  where  w  <•>  w.,  then  the  correlation  is  near  unity,  as 


as  shown  by  the  peak  of  Figure  4,  for  x  <  w^.  If  the  projections  are 
as  in  Figure  6b,  then  the  correlation  is  as  shown  in  Figure  4  between 
A  and  B  The  correlation  is  then  proportional  to  the  degree  of  overlap 
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of  strips  of  width  along  the  edges  of  the  two  projections.  These 
overlapping  strips,  as  x  increases,  move  like  the  slides  of  a  trom¬ 
bone.  If  finally  the  projections  are  as  in  Figure  6c  ,  then  turbulence 
in  the  range  of  strip  overlap  will  deflect  light  into  one  detector  when 
it  deflects  light  out  of  the  other,  which  gives  rise  to  the  anti¬ 
correlation  shown  in  Figure  4  for  x'*'  Zx 


'4 


Figure  7 

When  the  line  of  centers  of  the  detectors  is  quite  non-parallel 
with  the  source  sides,  theh  the  projections  of  interest  are  as  shown 
m  Figure  7.  There  is  then  no  extensive  overlap  of  the  edge  strip 
which  yields  the  correlation.  The  trombone  effect  is  then  absent,  and 
correlatior  varies  with  detector  separation  as  shown  in  Figure  5. 

From  the  above  results  for  a  rectangular  source,  it  is  possible 
to  closely  estimate  the  correlation  as  a  function  of  detector  separa¬ 
tion  when  the  source  is  the  sun.  In  this  case,  when  the  projections 
of  the  sun  are  as  in  Figure  8a,  with  w  <  Wj  ,  then  <r  n  varies  with 
x  essentially  as  in  Figure  5.  Eut  because  of  the  circular  nature  of 
the  sun,  cr^  will  fall  off  somewhat  less  rapidly  for  w  >  w^  than  in 
Figure  5.  When  the  projections  of  the  sun  on  the  turbulent  layer  are 
as  in  Figure  8c,  the  correlation  will  be  essentially  zero.  When  the 
solar  projections  are  as  irvFigure  8c  ,  there  will  be  a  slight  anti- 


.  *  —  —  aid  in  <JC  ,  out  not  so  promi¬ 

nent  as  if  the  source  were  square  with  two  sides  parallel  to  the 
line  of  centers  of  the  detectors. 

The  above  discussion  shows  that  the  anticorrelation  which 
occurs  for  a  rectangular  Source  with  specially  placed  detectors,  or 
to  a  lesser  degree  with  the  sun,  depends  on  the  height  of  the  turbu¬ 
lence.  In  fact,  the  observation  of  the  detector  separation  at  which 
the  anticorrelation  occurs  can  be  used  in  principle  to  measure  the 
distribution  in  altitude  of  turbulence.  If  a  square  light  source  were 
carxiec  on  a  high  altitude  balloon  and  observed  from  the  ground,  the 
method  might  be  useful  if  the  distribution  of  turbulence  were  suf¬ 
ficiently  simple,  that  is,  |f  it  consisted  principally  of  only  two  or 
three  well  defined  layers. 

When  light  from  the  sun  falls  on  two  nearby  small  detectors, 
the  correlation  between  fluctuations  in  their  outputs,  as  a  function  of 
their  separation,  will  be  a  single  hump  of  half  width  of  the  order  of 
STXT  where  L  is  the  mean  height  of  the  turbulence.  Equipment  of  ver 
great  accuracy  could  in  principle  show  the  details  of  the  edge  of  this 
hump,  associated  with  the  circular  nature  of  the  solar  source,  and 
could  in  principle  show  the  anticorrelation  when  the  angle  subtended 
at  a  point  in  the  turbulence  by  a  line  between  the  sources  becomes 
equal  to  the  angle  subtended  at  a  point  in  the  turbulence  by  a  solar 
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®*  .^■’^e  Absolute  Value  of  Fluctuations 

The  preceding  sections,  except  for  a  factor,  solve  the  problem 
of  the  correlation  of  fluctuations  in  detector  currents  produced  by 
the  light  from  an  extended  source.  The  results  have  been  expressed 
in  terms  of  a  correlation  coefficient^  (t),  which  is  unity  when  the  de¬ 
tector  separation  b  is  zero.  To  be  complete,  these  results  must 
all  be  multiplied  by  P,  the  rms  value  of  the  fluctuations  per  unit  area 
of  detector,  when  the  source  is  so  small  that  p( b)~l  where  b  is  equal 
to  a  source  diameter.  Of  course,  they  must  also  be  multiplied  by 
detector  area  B  to  give  proper  absolute  values. 

The  value  of  P  can  be  computed  from  formula'  (4.  18)  of  Ref.  1, 
namely. 


dS  = 


ck4A  2 sin2 
o 


4z 


—  (k-km)dV 


(43) 


where  dS  is  the  mean  scattered  intensity  at  distance  from  a  volume 
dV  containing  turbulent  eddies  of  air.  c  is  the  velocity  of  light,  £  the 
incident  propagation  vector,  Aq  the  electric  field  in  the  incident 
light  beam  prior  to  scattering  by  turbulence,  while  inis  a  unit  vector 
m  the  direction  from  dV  to  the  detector.  %  is  the  angle  between  k 
and  m.  0n  is  the  Fourier  space  component  of  the  unnormalized  corre¬ 
lation  B(r)  of  the  index  of  refraction  n  =  1  +  within  the  turbulence. 
More  completely  stated 


B(r)  =  n^Pjn^r'+r) 


(44) 


where  the  bar  indicates  a  statistical  average,  and  nj  is  the  fluctuation 
in  the  index  of  refraction.  Also,  by  definition  of 


B(r) 


r 


j4’(k)  die 
n 


(45) 


where  k  «  Z-x/A  is  the  length  of  the  wave  vector  of  the  spatial  Fourier 
component  of  iij  with  the  wavelength  A.  Tatarski  shows  that 


(46) 
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0(k)  =  £(8/3)  Sin  (vr/3)  2  -  l  i/3 

a  TT  —  c„  «  ' 


4ir 

2T^0[[moial  OJthc  spoctrum  of  tarbulence  u  ueed  • 
and  when  Cn  expresses  the  strength  of  the  turbulence.  C  ,  propo-_ 

IT!  '°  r  ValUe  ni<r)'  m"St  b°  — •»«**»—  from  experiment 
“  °rde-  °  °btain  *»  «  «*.!  atmosphere.  Evaluating  the 

numerical  coefficients  in  (46)  yields 


p(x)  =  0,  033  C  2  k"11/3 
n  » 


*o  <  K  <  *i 


(47) 


iocated  i  7;  °XpreSS°S  'here  CXiS‘  eddUs  °f  diameter  w 
located  m  the  ft  eld  of  turbulence  end  whose  diffraction  half  angle  tfw 

“  :r  r angie  eub,ended  by  ^ « *•  — .  d, 

iron  is  usually  satisfied  for  turbulence  in  a  clear  atmosphere, 
under  ih"  ^  ^  ^  °btained  from  data  on  s‘ar  twinkling 

rrrr* turbuience  is  u“  *• ™ * 

h,;‘f  .  °  *“•  Wth  a  mea"  »"***  of  10  km.  A  consistent  value 

ft  s  a  so  been  obtained  from  radar  scattering  experiments.  The 
result  is  „ 


0> 


C  -  7  .  lO^cm-1/3 


(48) 


The  value  of  Cr  is  related  to  t 
the  equation 


emperature  fluctuations  in  the  air  through 


Cn  =  79  •  ,10"6-^ 


pc/: 


who: 


p  is  pressure,  T  is  absolute  temperature  and  C.p  is  defined 


[T(r+*i)  -  T^)J  2  =  CTr2/3 


(49) 

by 

(50) 


The  quantity  on  the 'left  offSOJi,  called  DT(r,,  the  structure  function 
of  the  temperature.  A  simple  algebraic  expansion  show. 


>  »  t.  ,  ^ 
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DT(r)  *  2[t^..  Tjtr+rj)  T(r)j 

where  Tj  -  T  -  T  is  the  temperature  fluctuation,  and  D  is  twice 
the  squared  variance  minus  the  correlation  of  T,.  The  value  of  C 
consistent  with  (48)  is  8  •  10“ 3.  T 

The  angle  Xia  (  t)  is  given  by 
cos  %.  =  sin  0  cos  fi 

•  ^ 

where  0  is  the  angle  between  the  propagation  vector  and  a  line  con¬ 
necting  the  detector  to  the  dV  and  fl  is  the  azimuth  of  A  .  Scattering 
angles  are  of  the  order  of  w/z,  .  where  w  is  the  scale  of  turbulence 
and  zt  its  altitude,  in  a  clear  atmosphere,  these  are  so  small  that 
we  put  sin%=  i.  By  simple  geometry 


£  -  kin  a  2k  sin  0/2  -  kq/z 


where  q  is  distance  of  dV  from  the  line  connecting  the  source  point  and 

the  detector. Also,  dV  .  hdq  ,  where  h  is  the  thickness  of  the  turbulent 
layer.  Substituting  in  (1)  yields 


where 


dS 


,4  i 

-y  0.033  G 


■  2/*r! 11/3  v, 

’n  (Kql  h  2Trq  dq 


So  s  A0  c/8lr 

is  the  Poynting  vector  for  unscattered  light  from  the  distant  source. 
Strictly,  q  should  be  taken  as  a  function  of  altitude;  but  for  the  present 
we  evaluate  it  at  the  mean  height  of  turbulence,  namely,  zy  *  10  km. 

The  quantity  (54)  is  to  be  integrated  over  the  range  q  <  q  < 
where  qQ  =  (Azr)  IS  the,  diameter  of  the  smallest  eddy  which  contri¬ 
butes  strongly  to  fluctuations  at  the  detector  in  the  diffraction  zone;  and 

light.  We  have 
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By  (54)  and  (55) 

(5/So)  (2*)7/3  0.  033  Cn2h  (z,)5/6/^'7/6 


(55) 


(56) 


We  take 


f 


Si  c-  7  *  10  ^cm  ^  ,  h  a  20  km,  z  a  10  km  , 

r 

A=  5  •  10-5  cm 

.  ! 

Then  (14)  becomes 

\ 

’  \ 

S/So  .  1. 4  «  t  (57! 

Ty  4  v 

which  says  that  the  light  from  a  star  scattered  by  turbulence  into  a 
small  detector,  such  as  the  eye,  is  as  great  as  the  direct  light.  Since 
the  scattered  light  is  fluctuating  with  a  sero  mean,  this  corresponds 
to  strong  twinkling  of  a  star,  which  is  observed.  That  it  does  not 
appear  stronger  to  the  eye  than  it  does  may  be  caused  by  the  eye's  per¬ 
sistence  of  vision  (0~.  \  sec)  and  the  fact  that  the  mean  frequency 
of  the  fluctuations  of  the  scattered  light  is  about  100  cps  except  when 
wind  speeds  are  unusually  low.  The  excess  of  S/Sq  over  1  cannot  be  taken 
seriously  and  even  contradicts  conservation  of  energy.  The  above 
approximate  calculation  is  correct  only  to  a  factor  of  about  2.  In 

what  follows,  we  take  S/S  -  1. 

9  o 

3 

The  variance  <r  of  turbulence  produced  by  fluctuations  from  the 
sun  can  now  be  c  omputed  from  (37).  We  take 
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Do 


As  already  mentioned  above,  a  detailed  correlation  coeffi¬ 
cient  for  fluctuations  at  the  ground  of  light  received  from  a  point 
source  is  given  by  Tatarski  in  Reference  1.  He  obtains  a  theoretical 
correlation  coefficient  (p.  154)  and  also  an  experimental  one  (p.  213), 
which  if  reproduced  in  the  upper  portion  of  Figure  9. 


An  approximation  to  the  experimental  curve  (solid  line)  is 
given  in  the  form  of  Equation  (21)  ^ash  line)  with  c^  =  1.  33,  =  0.  33, 

ui  “  0.  956  YaL'  and  U2  a  1 .  9 1 2 tA L  .  These  numbers  assure  that 
the  Condition  22  is  fulfilled.  It  should  be  noted  that  the  experimental 
curve  does  not  satisfy  the  Condition  22,  which,  as  discussed  above, 
is  necessary  when  the  turbulent  atmosphere  is  non-absorbing  and  a 
sufficiently  weak  scatterer.  It  is  clear,  however,  that  the  solid  line 
could  be  shifted  considerably  within  the  indicated  experimental  error 
so  that  the  Condition  22  would  be  satisfied. 


The  central  curve  in  Figure  9  shows  the  correlation  coefficient 
for  the  intensity  fluctuations  in  two  point  detectors  as  a  function  of 
their  separation  distance.  This  is  obtained  by  assuming  that  the  Tie 
of  separation  between  the  two  detectors  lief  along  the  diagonal  of  the 
assumed  square  source.  It  is  seen  that  the  half  width  of  the  intensity 
correlation  is  approximately  l^L  ^7  cm. 

The 'trombone  effect"  is  demonstrated  in  the  lower  portion  of 
Figure  9.  The  change  to  a  logarithmic  scale  should  be  recognized. 
This  curve  shows  the  correlation  over  large  distances  when  the 
detectors  lie  on  a  line  corresponding  to  a  side  of  the  square  source, 
ibis  curve  demonstrates  clearly  the  edge  effects  as  discussed  above. 
The  anticorrelation  dip  at  b/Y/TL  ~2  x  103  is  obtained  for  a  source  of 
solar  dimensions  and  occurs  at  approximately  the  distance  of  the  pro¬ 
jection  of  the  sun  onto  the  earth's  surface  through  the  turbulent  layer, 
""140  xru  The  width  of  the  central  peak  and  the  width  of  the  anticorre¬ 
lation  dip,  are  approximately  equal  to  VJL, although  at  first  glance  this 
is  not  apparent  because  of  the  change  of  scale. 


Empirical  correlation  coefficient  of  tho  fluctuations  of  light 
intensity  (Ref.  1).  ) 

Approximate  correlation  coefficient. 


Correlation  coefficient  for  extended  source  with  point  detectors  along  diagonal 
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APPENDIX  E 

Experimental  Data  on  Scintillation  of  Sunlight 
By  R.  Hopfield 

We  have  experimentally  studied  the  fluctuations  of  sunlight 
caused  by  atmospheric  scintillation.  Preliminary  results  have 
been  obtained  concerning 


■I 


1.  Absolute  noise  levels  due  to  scintillation. 

2.  Correlation  of  signals  to  separate  detectors. 

3.  Power  spectra  of  the  noise. 

Some  of  these  Sadies  were  made  during  inversion  conditions,  a  not 
uncommon  situation  in  the  Los  Angeles  area.  The  effects  of  these 
conditions  on  the  experimental  results  are  not  fully  understood  at 
the  present  time.  Phenomena  not  related  to  the  refractive  effect 
of  turbulence  may  affect  the  results  of  experiments  made  under  these 
conditions,  i.  e. ,  effects  such  as  structured  haze  blowing  by,  vari¬ 
ations  in  the  water  vapor  in  the  path,  etc.  All  quantitative  data  given 
were  taken  with  IRC  SI 020E8- solar  cells  (area  2  cm2)  individually 
biased  at  10  ohms.  Multiple  cell  configurations  were  formed  by 
series  connection  of  the  biased  units. 

Absolute  Noise  Levels 


The  absolute  noise  levels  are  given  in  Table  I.  The  data  were 
taken  in  the  Santa  Monica  area.  Noise  level  varies  considerably 
from  day  to  day  and  hour»to  hour.  No  attempt  has  been  made  to  arrive 
at  an  average;  the  figures  given  in  Table  I  are  only  representative. 

The  "closed"  configuration  was  a  5  x  5  array  of  cells  with  overall 
dimensions  7x10  cm;  the  "open"  configuration  also  had  five  cells  to 

a  side,  but  the  cells  were  mounted  with  a  three-inch  spacing  between 
centers. 

f  ■  ,  '  V  jv 
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2.  Signal  Correlation 

Figure  1  shows  actual  unselected  scope  traces  generated  by 
wo  cells  at  various  separations.  The  sense  of  the  lower  tra-e  is 
opposite  to  that  of  the  upper  trace,  to  allow  the  use  of  symmetry  in 
evaluating  correlation.  It  is  obvious  that  correlation  exists,  and 
a  it  decreases  with  increasing  detector  separation,  but  we  have 
not  yet  succeeded  in  getting  any  trustworthy  numbers  in  this  area. 
The  data  shown  were  taken  during  inversion  conditions. 

3'  ~2^g£  SpV  fa  of  Scintillation  Noise 

M  „  ,Tihr:  Wre  CViUUated  With  4  «”*«  Panoramic  Analyser. 

°  !«  "  A"  ThC  mea6Ured  'quare  bandwidth  of  the  analyzer  was 
“  CyC,a*‘  FigUre  2  6hows  'Pactral  data  for  the  output  of  1 
on  .  cells;  integrated  noise  in  each  case  corresponds  within  30 

percent  to  the  total  noise  at  the  time  measured  on  a  Ballantyne  true 

rm.  meter.  By  inspection  of  Figure  2  it  appears  that  roughly 

J  **  t0*al  “^iUation  «**>  P°wer  is  in  frequencies 

°7  CPS'  thC  b«»g  o  little  more  than  a  factor  of  10 

per  100  cps. 


TABLE  l 


Absolute  Scintillation  Noise  Level 


Close  Array  -  No  Inversion 


No.  of  Cells 


CellAraa  DC  Voltage  rms  Noise 

i  cm  )  /  t 


Noise  Powc 
M  watts 


(volts) 


Open  Arra 


No  Inversion 


No.  of  Cells 


rms  Noise 
(volts) 


Noise  Power 
M  watts 


20  m  o  c/cm 


Power  Spectra  of  Scintillation 


r  inversion  conditions,  where  the  origin  of 
ib  not  fully  understood. ) 


